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We consider a system of particles in an array of microscopic traps, coupled to each other via elec- 
trostatic interaction, and pushed by an external state-dependent force. We show how to implement 
a two-qubit quantum gate between two such particles with a high fidelity. 
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I. INTRODUCTION 



The possibilities offered by quantum mechanical sys- 
tems for efficient information processing have stimulated 
in recent years the rise of an entirely new field of re- 
search [Q . Quantum protocols for secure communication 
over long distances have been devised and demonstrated. 
Quantum algorithms for efficient solution of problems be- 
lieved to be intractable on classical computers have been 
developed. However, while quantum communication is 
already approaching the stage of real- world applications, 
quantum computation remains still at a less advanced 
level, as far as physical implementation is concerned. Dif- 
ferent systems are being proposed as candidates for this 
purpose [^| , but nobody can yet tell what will turn out to 
be a viable solution. Indeed, in a few cases quantum com- 
putation building blocks - single- and two-qubit opera- 
tions - have been already demonstrated experimentally. 
In principle, these ingredients are universal - they are suf- 
ficient to build an arbitrary unitary transformation over 
N qubits (i.e., any quantum computation). But in order 
to perform useful computations in a real environment in- 
ducing decoherence, fault tolerance is also required. This 
implies e.g. nested redundant coding for real-time error 
correction |}|, and requires an error probability for el- 
ementary operations below a certain threshold (of the 
order of 1CP 4 ). Hence the need for new proposals, al- 
lowing for handling a bigger number of qubits at a lower 
decoherence rate and with faster and more reliable gate 
operations - in a word, enabling scalability of the system. 

We propose to use quantum optical systems in periodic 
microscopic potentials. This is meant to combine the 
good isolation and precise control by laser fields, achiev- 
able in quantum optics, with the ability - usually associ- 
ated with semiconductor technology - of manufacturing 
periodic structures to generate modulated fields on a mi- 
croscopic scale. The general concept of our proposal is 
to encode the logical states of each qubit into two in- 
ternal states of a particle (neutral atom or ion). Single- 
qubit operations are obtained as Rabi rotations by apply- 
ing resonant laser fields. Two-qubit gates are performed 
by inducing a state-dependent interaction over a certain 
time, making the particles acquire a conditional phase 
shift depending on their logical states. These, however, 



in a real situation are coupled to other external degrees 
of freedom. This can lead to different kinds of imper- 
fections. On one hand, the external state after gate op- 
eration may not be exactly the same as before. On the 
other hand, the conditional phase shift will also depend 
on the external state: if this is mixed, only an imprecise 
phase determination will be possible. These facts affect 
the gate fidelity, which is defined by comparing the de- 
sired effect of the gate with the actual evolution that can 
be obtained in the laboratory. We already proposed sev- 
eral schemes, based on different interactions - collisional 
interactions between neutral atoms in optical lattices Q| 
and magnetic microtraps |j| , or dipole-dipole interactions 
between Rydberg-excited atoms Hj. Here we deal with 
electrostatic interaction between ions [0, || in arrays of 
microscopic traps Q. 

In this paper we describe the conditional dynamics 
for two charged particles, trapped in separate harmonic 
wells, interacting via electrostatic repulsion and under 
the influence of an external state-dependent force, which 
can be generated e.g. by an off-resonant laser standing 
wave |8| . The goal is to implement a phase gate between 
the two qubits (Sec. II), i.e. to transform their initial 
state by inducing a certain phase onto each of its compo- 
nents. The ideal transformed state so defined has to be 
compared to the one that can be obtained by means of 
a realistic Hamiltonian, coupling the par ticles' internal 
and external degrees of freedom (Sec. III). To this aim, 
we consider first a one-dimensional classical model for 
the motion (Sec. IV). We solve the equations of motion 



for each combination of logical states separately, and de- 
fine in each case a two-particle phase as the integral over 
time of the interaction energy. These phases define the 
evolved internal state, to be eventually compared with 
the ideal state we aim to obtain. The fidelity of gate op- 
eration is evaluated, in this classical model, as the overlap 
between the real and the ideal state, averaged over dif- 
ferent starting conditions according to a thermal proba- 
bility distribution for the initial oscillation energies. The 
outcome is a series expansion for the temperature depen- 
dence of the fidelity, of which we give explicitly the first 
terms. The full three-dimensional quantum-mechanical 
calculation follows the same path (Sec. [v]), except for a 
couple of points. First, the two-particle phases are now 
given directly by the Schrodinger equation and not de- 
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fined ad hoc as before. Second, the fidelity is evaluated 
by tracing out the external variables of a density ma- 
trix representing a mixed thermal quantum state. We 
calculate perturbative corrections arising from multipolc 
terms in the Coulomb potential, and show how to sup- 
press lowest-order corrections to the fidelity by means of 
an intermediate tt rotation on the qubits, thus achieving 
an improvement by several orders of magnitude. 



II. A QUANTUM PHASE GATE 

We want to implement quantum logic between parti- 
cles stored in an array of microscopic traps. The qubits' 
logical states |0) and |1) are encoded into particles' in- 
ternal states. One basic building block towards multi- 
qubit entanglement operations is the phase gate between 
two qubits - a transformation which rotates by a certain 
phase just one component of logical states: 
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When i9 = 7r, this is equivalent - up to single qubit rota- 
tions - to a Controlled-NOT gate. Ideally, this would be 
accomplished by means of a state-dependent interaction 
of the form 



fftot = A£7(t)|l)i(l| ® |l) a (l|, 
acting over a time r such that 

AE(t')dt' = i9. 



(2) 



(3) 



However, it is not straightforward to realize in practice 
an interaction between two particles which couples only 
their internal states - other degrees of freedom, for in- 
stance the motional ones, are likely to be affected. There- 
fore our goal is to approximate the ideal transformation 
Eq. (^) by means of a conditional dynamics for two parti- 
cles, making them acquire the phase •& if and only if they 
are both in the internal state |1), and leaving eventually 
the external degrees of freedom practically unaffected. 
This is described a Hamiltonian of the form 



JJ(t,xi,x 2 )= J2 # a/3 (*.xi,x 2 ) Hl< 

a, ,3=0 



(4) 



where Xj denotes the external degrees of freedom of par- 
ticle j, and the explicit time dependence indicates that 
we can switch on and off a suitable interaction in order 
to obtain the desired effect. To evaluate the performance 
of our scheme, we have to compare the case of an ideal 
gate, as given by Eq. p), with the gate that can be actu- 
ally realized by the physical process described by Eq. (El). 



The figure of merit is the minimum fidelity F, given by 

(5) 



F = min tr ext (xVlx') 
x 



where tr ext denotes the trace over the external degrees 
of freedom, \x) = c a/3\ a )i 1/3)2 is a generic two-ion 
internal state, |x') is the state obtained from \\) via the 
transformation Eq. (|l|), and a' is the total density matrix, 
including external degrees of freedom, after the evolution 
dictated by the Hamiltonian Eq. (Q), starting from an 
initial state 



a = pi (t ) <8>p 2 (*o) ® |x)(xl) 



(6) 



where Pj(io) is the external state of particle j at the 
initial time to. Ideally, to achieve the optimal fidelity 
F = 1, we need that the external degrees of freedom 
factorize after the gate operation, and that the evolution 
operator 



U(t,t ) 



(7) 



has the only effect to induce a two-particle phase <p a P , 
depending on the internal state of both ions, plus single- 
particle phases, due to the kinetic energy associated with 
the trap displacement. The latter can be undone by 
means of single-qubit rotations (see App. Bl), leaving 
us with the gate phase 



(8) 



In a real situation, the starting point will be rather 
a mixed state corresponding to a thermal distribution 
over the external energy eigenstates. In other words, at 
nonzero temperatures there will be a finite probability 
that each particle starts in an excited motional state, 
leading in general to different phases, which cannot be 
experimentally controlled and easily undone by single- 
qubit rotations. Therefore the fidelity is expected to de- 
crease with temperature, as we are going to show quanti- 
tatively in the next Sections, both in a classical model for 
the particles' motion and in a fully quantum framework. 



III. CONDITIONAL DYNAMICS 

We consider N ions, trapped at positions denoted by 
(c-numbers) Fj (1 < i < N). For simplicity, we take 
the trapping potentials for all ions to be harmonic, with 
the same frequency u) along every spatial direction. Our 
results can be straightforwardly generalized to inhomo- 
geneous trap arrays with anisotropic confinement. More- 
over, each ion is assumed to be subject to a time- varying 
force Fi(t), depending on its internal state on € {0, 1} as 
in Eq. (||). The Hamiltonian is 



N 



N 



(9) 
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where 



H t = ^ + Jmw 2 (r 1 -f t ) 2 -F 1 (t) -r,, (10a) 
2m 2 



H, 



(17 
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In the following, we will focus on two-particle dynamics. 
We assume the external force to have the same strength 
on both ions, i.e. to depend only on the internal state of 
each particle: 



F,(t) ®F«(t). 



(11) 



We can rewrite the Hamiltonian in Eq. (0) as 



ff aiaa (t,xi,x 2 ) 



E 



i=l 



pi muS 

2m 2 1 



x"*W] 2 (12) 



F Q -(i) 2 
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Atteq Id + x 2 - xi| 



where we have defined xi = ri + d/2, x 2 = r 2 — d/2, 
x a (t) = F a (t) / (mu> 2 ) , and d is the equilibrium inter- 
particle separation. We assume the external force to act 
along its direction, and we choose the simple state depen- 
dence F a (t) = (ahu>F(t)/au, 0, 0). Here a w = ^h/miv 
is the quantum harmonic oscillator ground-state width. 
Hence the adimensional quantity J-(t) represents the dis- 
placement, in units of a Ul induced by the force on the trap 
minimum for ion i, if it is in internal state |l)j: indeed, 
we have x Q (t) = aa u (!F(t), 0, 0) . With the above choice 
for the state dependence of the force, the last term in the 
first row at right-hand side of Eq. (|l^) will not contribute 
to the gate phase Eq. (||) , since the corresponding terms 
cancel each others in the sum z~2 a p(~l) a+l3 <p al3 ■ Since 
the interaction only depends on the distance between the 
particles, it will affect only the relative motion. Therefore 
we can study the problem in the coordinate system where 
the relative motion is decoupled from the center-of-mass 
degrees of freedom. The Hamiltonian can be rewritten 
(see App. lO) as H(t) = H R (t) + H r (t), where 



H R (t) = Hi - F(t) • R 



ri + r 2 



H r (t) = H° r -f{t) -(r + dj+ffi. 



(13a) 
(13b) 



Here, H R and contain three-dimensional harmonic 
potentials, and describe also nonadiabatic effects arising 
when lot ~ 1. In particular, H® incorporates terms aris- 
ing from the interaction up to the order (a^/d) 2 . Hi 
entails the higher-order multipole contributions. 



IV. A CLASSICAL MODEL 

We first treat the ions' motion classically, i.e. we re- 
gard them as point particles following well-defined trajec- 
tories dictated by the state-dependent trapping potential 
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FIG. 1: Gate operation scheme in a classical picture. Above: 
selective trap displacement for the two different internal 
states. Below: schematics of the four different Hamiltonians 
for each combination of internal states. 



and by the repulsive electrostatic force. Without loss of 
generality, we can take the x axis parallel to d. We will 
study the one-dimensional problem of the motion along 
that direction, denoting by italic letters the first Carte- 
sian component of the vectors defined in the previous 
Section. The initial state of the system is described by 
the internal quantum state 



IXci(to)> 



E C "/3l a )ll 
a,P=0 



(14) 



and by the external classical trajectories it) of the 
two ions, dictated by the Hamiltonian Eq. ([|). Here, by 
a^aj (f3 f± a 2 ) we mean the internal state of the first 
(second) particle. Indeed, to find the trajectories for all 
values of a and (3, we have to solve four distinct classi- 
cal two-particle equations of motion, each describing the 
dynamics for one of the possible combinations of internal 
states, as depicted in Fig. (|l|). Once we have done that, 
we can evaluate the Coulomb interaction energy 
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We then define the evolved internal state as 

*3 



| Xcl (t)) = ^c^|a)|/3)e^S , 

a/3 



where 



a/3 
Vol 



(15) 



(16) 



(17) 



to 



Now we make the following assumptions: (i) the force 
acts slowly over the harmonic oscillator time scale 
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i.e. \T\ <§; oj; (ii) it induces a displacement of the or- 
der of the single-trap harmonic oscillator length a w ; (iii) 
the latter is much smaller than the distance between the 
traps, i.e. a w <C d; (iv) the amplitude of the intra-well 
oscillations (if any) is negligible with respect to the inter- 
well distance, i.e. Ei <C mLu 2 d 2 /2. The first three con- 
ditions can be fulfilled by construction, the last one re- 
quires in principle the motion to be cooled. Assumption 
(i) amounts to neglecting non-adiabatic terms in the tra- 
jectories (e.g., sloshing motion excited by the trap dis- 
placement). On the other hand, when the last three as- 
sumptions hold, we can consider, to a first approxima- 
tion, the intra-well motion to be basically unaffected by 
the higher-order multipole terms in the expansion of the 
Coulomb interaction, Eq. (fl5l). This approximation is 
not easy to check classically, since the exact trajectories 
cannot be computed analytically. We will test its valid- 
ity in second-order perturbation theory, in the con text of 
the quantum-mechanical treatment (see App. B 4). 



A. Starting conditions 

The initial classical motional state, at t — to, can be 
either the ground state, described by the initial condi- 
tions Xi(to) = ij(^o) = 0, or an excited state, described 
by oscillations of each ion inside its trap with an energy 
Ei (i = 1,2), i.e. by the initial conditions 



Si (to) = Aarf'(to) 



2E- 

cos [ui(t - U)} , (18a) 

mi . 



2Ei 



Xi(h) = -\l = — L sm[ui(t - ti)} 
w m 



(18b) 



(of course, the former is a particular case of the latter, 
for Ei = E 2 = 0). Here, w is a corrected trap frequency, 
taking into account up to quadratic terms in the Coulomb 
potential - i.e., up to n = 2 in Eq. ( |l5| ) - namely, Co = 
wyl + e/2, where 



(19) 



TTSomL} 2 d 3 



is essentially twice the ratio of the Coulomb energy 
q 2 /(4:Treod) and the energy of the second ion with respect 
to the first trap muj 2 d 2 /2. Under the approximations 



discussed above, we can write the trajectories x^ (t) as 



(20) 



at all times. The situation is depicted in Fig. |2| Note 
that we are treating classically the particle motion, but 
not the internal state: so the ions are allowed to be in 
a superposition of the available logical states, i.e. to os- 
cillate according so to speak to two different trapping 
potentials, as seen in Fig. || 
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FIG. 2: Gate operation dynamics for two classical particles 
oscillating with energies Ei. The state-selective trap dis- 



placements x ai (t) and the intra-well oscillations Ax t 
shown. 



B. Gate phases 



We can write the phases Eq. ( |i7|) as 



(f) are 



a/3 
Vol 



h<*0 



where 



is the ground-state contribution, and 



(21) 



is the correction due to motional excitations. To evalu- 
ate the various contributions explicitly under the above 
approximations, we now need only to specify the time de- 
pendence of the trap displacement. We choose the Gaus- 
sian form 



(22) 



Let us first consider the case where both particles are 
in their motional ground state. We insert Eq. ( |2p| ) for 
Ej = 0, through Eq. @, into Eq. @, and obtain 



-{a -(if 



i 2 tC0T 



OO 

n=0 



(23) 



In the evaluation of the ground-state phase <j>"f , Eq. J23| ) , 
the complete Coulomb potential Eq. ( |l5| ) has been taken 
into account. When we evaluate the corrections 8<t> a ^ 
instead, it is not possible to find a general expression valid 
at all orders n, which therefore have to be considered 
separately. We choose to as an integer multiple of the 
oscillation period 2ir/Co (so that the motional state is left 
unchanged after gate operation), and find 



A"/ 3 - 



= 3(a-/3)W-£W 



—= — - + (a — p) — 

V2£ d K 1 '\ d 



x j^r{Ei + E 2 + 2-s/ExE 2 cos[Q{t x - t 2 )]} 
+ o[(a /d) 3 }. (24) 

where it has been taken into account that lot 1. The 
two terms under square brackets in Eq. (^4|) come from 
terms in the Coulomb potential with n = 3 and n = 4 
in Eq. (|l5|), respectively. This means that no thermal 
correction is to be expected if only harmonic contribu- 
tions to the potential (i.e., with n < 2) are included. 
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Indeed, in this case the spurious interaction phases, due 
to the oscillations in the ions' positions, are averaged out 
when integrating on a time much larger than the oscilla- 
tion period. This explains intuitively why the phase does 
not depend on the motional state, in the approximation 
where only linear and quadratic terms in the Coulomb 
potential are tak en in to account, as it will be shown an- 
alytically in Sec . |VA| . Now, the classical analogue of the 
gate phase Eq. IM) can be written as 







a,8 



■Lii/2 [{£,a u /df 



where YA n (z) 
tion, and 



8 {a^/df 

e d + o[( aCj /d)% (25) 
Tl^kLi z k /k n is the polylogarithm func- 



(26) 



From Eqs. (E5h and fl26|), we see that in our one- 
dimensional classical model, up to first order in powers of 

/d, the gate phase is insensitive to motio nal excitations 
inside each trap. Later on (see Sec. VE ), we will find 
that the very same expression for the gate phase can be 
obtained, under the same approximations, with the full 
three-dimensional quantum formalism. Moreover, note 
that assumption (ii) applied to Eq. ( |22] ) means £ 1, 
whence (i) implies lot 1 . It follows that, if we want to 
obtain O \ = ir, it has to be e <C 1. This means that the 
confinement has to be strong with respect to Coulomb in- 
teraction over the inter-well separation, which is in turn 
consistent with assumption (iii). 

We will now consider a more general initial condition 
than the ones discussed so far, namely a thermal state, 
described by a probability distribution over the energies 
Ei and the oscillation phases u)t\. Assuming the energy 
distribution characteristic of a canonical ensemble and a 
uniform probability distribution for ujt\ and U)t2, we can 
compute the thermal averaged phase 



2n/u 



-0,1 <^ r 



+ 



\ d 



dt\dt2 
1 a- (3 



dE 1 dE 2 



2y/2 



d | ag, 
a& d 
6k B T 



n a P o k R T 



e | Qk B T 
<W + o[{acj/df 



(27) 



The mean gate phase turns out to be 



Indeed, we will see that the very same structure for the 
corrections to the lowest-order phase is obtained with the 
full quantum-mechanical calculation. 



C. Gate fidelity 

In order to obtain the desired phase gate Eq. (|l|), we 
require that ((i? c i)) = i". So the reference state, represent- 
ing the ideal evolution, is chosen as 



l^cl) = 



(29) 



a/3 



The real evolved state, Eq. (|T^), can be written as 
IXcl(t)) = |Xci) + I^Xd), whereby 



a/3 



(30) 



In our classical model, we are treating our particles' ex- 
ternal degrees of freedom classically. Therefore, in the 
evaluation of the fidelity Eq. (||), instead of tracing over 
the motional eigenstates we should average over the pos- 
sible classical trajectories. Thus 



F c x = rnm^( X ' cl \Xci(t))(Xci(t)\x' cl 



(30) 
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|l + <Xcl(*)|<5Xcl>r 
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a, f3=0 
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0\ e 



o[(a Q /d) 5 } 



; 2 V d) 



- 2 



(31) 



as is discussed in detail in App. [A 2\ . Finally, let us con- 
sider what would come out if we were able to suppress 
the cubic anharmonic correction from the Coulomb po- 



tential, i.e. to put k = 0. We will show later (Sect. VDQ 
how this can be done in practice - here we would like to 
give a classical estimate F c 'j of the im prov ed gate fidelity. 
The calculation is performed in App. A 2 as well, and the 
result is 



Ki(T) = 1 - 



W cX k B T 



(^Y + o[(a Q /d)% (32) 



This shows that, by suppressing one order of anharmonic 
corrections, one obtains an improvement by two orders 
in a^/d (several orders of magnitude) in the fidelity, as 
is shown in Fig. ^. 



a,0 



6k B T 



(28) 



i(ao/df 



V. QUANTUM TREATMENT 

We want to describe quantum mechanically the three- 
dimensional dynamics of the two particles. This means 
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FIG. 3: Fidelity F c \ (solid line) and improved fidelity 
(dashed line) in the classical model as a function of temper- 
ature T. Inset: detail of the departure from unity of the 
same quantities, on a logarithmic scale. We assumed to work 
with Ca + ions and chose the parameters uj — 2ir x 1MHz, 
d = 20/im. 




FIG. 4: Gate operation scheme in the quantum regime. Trap 
parameters are defined like in the classical case; the harmonic- 
oscillator ground-state width is also shown. 



that, unlike in the previous Section, their motional state 
is given by a wavefunction (see Fig- which evolves ac- 
cording to the Hamiltonian Eq. (jj). To better under- 
stand its structure, it is useful to write 

H a ?(t, xi , x 2 ) = ITf (t, x Q (t), x' 3 (f )) + Hf (t, xi , x 2 ) . 

(33) 

Here the Xj are the ion position operators, and the c- 
numbers x™' (t) denote the centers of the initial motional 
wavefunctions as determined by the trap (see Sect. III). 
To second order in the expansion Eq. ( |l5| ) , the first term 
at the right-hand side of Eq. ( p3| ) gives rise to the same 
contribution to the phase already calculated in lowest or- 
der in the classical model, namely 9 c i given in Eq. (p6|). 
Corrections to this phase are due, as before, to: (a) ther- 
mal excitations; (b) higher-order (multipole) terms in the 
expansion of the Coulomb potential; (c) nonadiabaticity. 
Motional effects of the kinds (a) and (c) are accounted 
for by , while multipole corrections get a contribu- 
tion also from H a @ (t, x Q (t) , x' 3 (t)) . In order to minimize 
such corrections, we choose to operate in the adiabatic 
regime, given by condition (i) in Sect. IV, i.e. we assume 
lot 1. We study the dynamics in the center-of-mass 



and relat ive coordinate system, as given by Eqs. (13a) 
and (13t). In both coordinate systems, the motion along 
different axes decouples: the transverse directions con- 
tribute just an overall phase, whereas the relevant state- 
dependent dynamics takes place along the x axis. Since 
we assumed d S> a u , it follows that H\ - containing only 
terms of o[(a u /(i) 3 ] - can be treated as a small perturba- 
tion in H . We will first neglect it and solve exactly the 
three-dimensional Schrodinger equation, and then take 
it into account perturbatively, eventually checking our 
results with a numerical simulation. 



A. Unperturbed forced oscillator 



When we put Hi 
solution to Eq. (|B4 
is done in App. |B 2 
this approximation, 



= and take Tit) as in Eq. (E2|), the 
can be written explicitly |1C|]. This 
The gate phase turns out to be, in 



(-1)°+^ = 2£ 2 [$(w) - $(wvTT^)] , (34) 



a,/3 



where 



$(w) 
if (Mo) 



as A(s, to)- 



ds 



(35a) 



i r 

= / dt' Tit')e lu{t - ta) . (35b) 

mhhJ Jt a 



Expli cit exp r ession s for (f> a ^ and for $(w) are given in 
Eqs. (|B12a )- (|B12(j ). In the limit wt > 1, we obtain 
<I>(w) w ~y / Tr/32LLiT. By expanding Eq. (|34| ) up to first 
order in e, we retrieve 9 w 9 C \ as given by Eq. (p6[). The 
phase 9 can be adjusted to the desired value 7r by tuning 
the displacement £ and/or the interaction time r. More- 
over, 9 is independent of the ions' motional state. This 
means that the phase remains the same even if we start 
with a mixed external state, described by a density ma- 
trix 



Prito) 



-H(t )/k B T 



(1-7) 6 



E 



7 '\n. 



= 1 ra; = l 



(36) 



corresponding to a thermal distribution at a temper- 
ature T. Here the canonical partition function Z = 
tr {e- H ^y kBT ), 7 = exp(-fi^/fc B T), {m}^,... , 6 = 
{n x ,ny,n z ,Nx,N Y ,Nz}, and the \N) X YZ i\n) x y z ) are 
the eigenstates of along each direction. To op- 

timize the gate fidelity (see Sect. II), we need basically 
one thing more - that the external degrees of freedom 
are not entangled with the internal ones after gate op- 
eration, i.e. that the final motional state does not de- 
pend on the logical states of the qubits. This indeed 
happens, under the adiabatic assumption (i) of Sect. IV. 
In fact, in this case, the overlap C^ ni ^(Mo) between the 
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initial and final spatial wavefunctions for a system start- 
ing in a motional eigenstate along all degrees of freedom, 
defined by Eq. (B14), is close to 1. To be more pre- 
cise, 0^ ni ^(t, to) formally depends on the motional state 
- however, if the adiabatic condition is satisfied and \t\, 
\to\ are large enough, 0' n ^i (t, tp) ten d s expo nentially to 
1, as it can be seen from Eqs. (B15a)-(B15b). 



B. Including higher-order terms 

We will now take into account the contribution of Hi . 
This does not affect the center-of-mass motion. There- 
fore from now on we will study only the relative motion. 
A little care is needed since, unlike in the most common 
time-dependent perturbation theory, in our case the un- 
perturbed Hamiltonian depends on time, while the per- 
turbation does not. The calculation is carried out in 
App. B3. In the adiabatic limit and for \t\, \to\ > t, first- 



order corrections simply amount to a relatively small ad- 
ditional phase shift: 



([/(Mo)) « (U (t,t ))e^ AC 



3 +A') 



Thus, to this order, 

ip al3 w 4> af} + A afj + A'. 
Hence we find $ = 9 + 58 + o[(a w /(i) 6 ], where 



(37) 



(38) 




t (S) 



x 10 



FIG. 5: Gate phase -ft/ir (circles), and projection of the initial 
motional ground state over the evolved one for ions in states 
1 0) i 1 1) 2 (stars), as induced by the external force J-{t) (dashed 
line). The results of the numerical calculation, performed 
with the parameters quoted in the text including corrections 
up to o[(a w ) ], are shown. Solid lines are analytical results in 
the harmonic approximation of Sect. VA. 



unity overlap of the final motional state on the initial 
one. Indeed, with t = —to = 150/is, even st artin g e.g. 
with the (10 4 ) th oscillator excited state, Eq. flB14| ) still 
predicts 0(t, to) > 1 — 1CP 10 . With these parameters, the 
perturbative estimate derived within the classical model 
in Sec. IV turns out to be C \ « 1.047T. 



V2£ 2 + 3(2n a 



n z , 



(39) 



and it is understood that we started from a pure state 
with n x ,y t z excitations along the various directions of the 
relative motion. As already anticipated, Eq. (^) has the 
same structure as the classical correction to the phase 
expressed by Eq. (|2^) , except for a different overall factor 
related to the dimensionality of the problem we are now 
considering. Indeed, the first term under square brackets 
comes from the higher multipoles (k = 3) in Hi, while 
the second one is due to the thermal excitations. 



C. Numerical computation 

In order to check the validity of the perturbative ex- 
pression Eq. ([39|) , we solved numerically the Schrodinger 
equation for the relative motion, taking into account cu- 
bic a nd qu artic inter action terms, explicitly given by 
Eqs. (A6a) and ( |A6b ) respectively. The calculation is 
described in App. |B E , and results are shown in Fig. |^. 
In particular we find that the cubic corrections cancel 
indeed each other, and that the quartic corrections have 
the same order of magnitude as predicted by Eq. (|39"|). 
From Eq. (|||) we obtain 8 = tt with id = 2tt x 1MHz, 
x 2 - xx = 20^m, £ = 0.7 and r = 41.1069/xs. These 
results were confirmed by the numerical computation up 
to 40 initial excitations in each direction, always giving 



D. Gate fidelity 

The gate phase $ cannot be measured directly, since 
the higher-order corrections arising from the Coulomb 
potential depend both on the internal and on the mo- 
tional state of each ion, and cannot be undone by means 
of single-qubit operations, unless the logical state is mea- 
sured. However, the corrections of order k have a simple 
internal -state dependence of the kind (a — (3) k , as shown 
by Eq. QB22a ). This implies that it is possible to obtain 
a cancellation of the odd-order corrections, by applying 
a 7T-pulse R = |0)(1| + |1)(0| to both qubits in the middle 
of gate operation. Indeed, if U is the evolution operator 
giving the dynamics described in the previous Sections, 
we find 



|0)|0) 
|0)|1) 
|1)|0) 

11)11) 



nu 



e^ °|l)|l) 
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,ltp 



1)10} 
|0)|1) 
|0)|0) 



— > e ^ 01 +^ 1Q )|o)|i) 

— > e^ 1O +^ O1 )|l)|0) 



(40) 



Here an adiabatic approximation is understood, accord- 
ing to which the final and initial motional state are iden- 
tical. We now define 



A8 = 



(41) 
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(as before, ((•)) denotes the thermal average), and the 
gate operator 



G = S{RUf 



(42) 



where 



e -2^ 2 *( w ) + | 1 \/ 1 | e -i[2« 2 *(a;VT+i)-<(5e»]_ 



If we choose the gate operation time r in such a way that 



7T = 4£ 2 [$(w) - $(wvTTe)] - 2((60)), (44) 
we obtain 



|0)|1) 
|1)|0) 
IDID 



-iA6> iG 



|0)|0) 
|0)|1) 
|1)|0) 

ll>ll>. 



(45) 



The global (thus irrelevant) phase 9 is given in App. B t . 
Note that the single-qubit rotation S is the same for both 
qubits, and therefore single-ion addressability is not re- 
quired. The fidelity of the gate operation is defined by 
comparing the ideal gate operation Eq. (|J) with the ac- 
tual dynamics obtained in our scheme, Eq. (^B]), at a 
given temperature T for the moti on i n all three dimen- 
sions. The result, derived in App. BE, is 



F(T) w 1 - 6 3 



(46) 



The fidelity turns out to be independent of r and £ which, 
subject to the conditions lot 3> 1 and £ ~ 1, can be freely 
chosen to obtain the desired gate phase. The dependence 
of the fidelity on the various parameters is the same as 
in the classical model discussed in the previous Section. 
As already anticipated in the previous Section, the inter- 
mediate 7r-pulse R allows us to get rid of the o[(a w /<i) 2 ] 
term, thus obtaining a much better gate performance. In- 
deed, the only difference between the corrected classical 
fidelity F' cl and the quantum fidelity F is the numerical 
pre-factor multiplying the temperature-dependent part, 
which is bigger in the latter case due to the inclusion of 
all the spatial degrees of freedom, whereas our classical 
model was just one-dimensional. Anyway, with the pa- 
rameters quoted above, at temperatures below 2 mK, cor- 
responding to an average number of harmonic-oscillator 
excitations n ~ 6, the fidelity turns out to be bigger than 
1 — 10 -6 . We can also evaluate how the fidelity scales 
when the gate is repeatedly applied, say g times. It is 
clear from Eq. ^ that in this case, apart again from an 
overall phase, 



|0)|0) 
|0)|1) 
|1)|0) 

11)11) 



G" 



e- l f Ae |0)|l) 
(-1) 9 |1)|1). 



(47) 



The excitation-dependent phase A9 is just multiplied by 
g. Thus, under the same approximations as above, the 
fidelity of the <?-fold gate operation is 

F^(T) = l-g 2 [l-F(T)], (48) 

i.e. it scales with the square of the number of gates. 

E. One-dimensional calculation for many ions 

We now assume to have N ions, trapped in a linear 
array of equally spaced traps, i.e. we take 



{x h y h Zj) =jd 



(49) 



in Eq. (|lOa|). Expanding the interaction Hamiltonian 
Hij, Eq. ( |10b| ), in powers of the i\ and rj, and neglecting 
terms of o[(a w /c?) 3 ], we find 



H 



N r 2 



E 



-±(xi - Xi) 2 + yf + zf 



to' 



2aJF{t)\l)i(l\{xi+Xi) 



N 

E 

Kj 



i-jf 



(50) 



where the various quantities are defined in App. B 7 



Eq. (|50|) describes a set of independent f orce d harmonic 
oscillators, like the ones we solve in App. B 2, plus a cou- 
pling term multiplied by e. If e <C 1, we can treat this 
term as a small perturbation, in the very same way we 
develope in App. B3. We take as initial state 



iV 



i**(*o)> = n i n ^i a * 



(51) 



where ctj denotes the internal state of the i th ion, and 
rii its motional state in the well corresponding to the i 
term of the first sum at the right-hand side of Eq. (50). 
We obtain 



I* 



TV 



J4>° 



I*jv(*o)), (52) 



where (calculating the two-particle phases <p aiaj pertur- 
batively) 



a 2 £ 2 <J>(w;) + aisfnLOT^Xi/ a u 

- [n x ,iLOi + (n y ,i + n z<i )io](t - t ), 



(53a) 



- ^ I ^N{h)\Uo{t,t')^\u {t',t )\^ N {t ))dt' 

t \ l 3\ 

2, 



\i-j\ 3 {l + er]i)(l + e Vj ) 



J 



;|3 



*cl, 



(53b) 



the last line following from e<l and Eq. ( B40 ). Again, 
the result to this order turns out to be independent of 
the motional state of anyone of the ions. In the case of 



two ions, Eq. (53b) gives back Eq. (26) 
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VI. CONCLUSIONS 

We analyzed in detail a recent proposal [|| for scal- 
able quantum computation with ions in an array of mi- 
crotraps. This scheme bears important advantages over 
the previous proposal JtJ, based on trapped ions as well. 
Indeed, in that case many ions, lying in a single trap min- 
imum, exchange information via the collective motional 
excitations; ground-state cooling is an absolute need, and 
any perturbation on each ion can affect the performance 
of the whole system. Here, instead, each ion is confined 
to a single minimum of a periodic microscopic potential, 
and interacts with other ions via Coulomb force. Un- 
der the conditions discussed in the text (adiabaticity of 
the trap displacement, strong confinement with respect 
to the distance between ions, intermediate symmetrizing 
7r-pulsc), the phase shift is insensitive, to a high accuracy, 
on the motional state of each ion inside each trap, and 
therefore the fidelity turns out to be practically indepen- 
dent from temperature. Moreover, trapping frequencies 



can be much higher than in the previous case, leading to 
much shorter gate operation times. As long as we take 
into account purely motional decoherence mechanisms, 
we find a fidelity bigger than 1 — 1CU 6 for a two-qubit 
phase gate operating on a time scale of a few tens of /is. 
Furthermore, with the improved scheme presented here, 
single-qubit addressability is not required for any of the 
various control operations. To sum up, the present pro- 
posal constitutes really a good candidate for a scalable 
implementation of a quantum computer. 
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APPENDIX A: CLASSICAL CALCULATION 

In this Appendix we give a detailed account of the cal- 
culations leading to the results obtained in the classical 
model for our two-qubit phase gate. 



1. Rewriting the Hamiltonian 

In this Section we give the explicit form of the vari- 
ous terms in Eqs. (13a) and (13t), describing the Hamil- 
tonian Eq. ([)]) for two ions, in the center-of-mass and 
relative-motion coordinate systems. Indeed the Hamilo- 
nian, with the number of ions N = 2, may be rewritten 



as H = Hr + H r , where 

p2 X 



H 



R 



2M 2 
„2 1 



Mu/R 2 - F(t) • (R + R ) (Ala) 



H r = y + |/^ 2 (r + d-d ) 2 
-f(t).(r + d)- 



A 



and 



R 



(A, Y, Z) 



ri + r 2 



Rn 



(x,y,z) = r 2 -r x - d, 



Rn 



ri + r 2 



Pi +P2, 



do 



r 2 -ri 



P2 - Pi 



(Alb) 

(A2a) 
(A2b) 

(do, 0,0), 

(A2c) 
(A2d) 



P(t)=Fi(t)+F 2 (t), f(t) 



F 2 (t)-Fi(i) 



M=2m, M= "2 ' 



A : 



4ire 



(A2e) 
(A2f) 



In the above, d = [d, 0,0) is the equilibrium separation 
between the two particles in the absence of the pushing 
force. Due to the repulsive Coulomb interaction, |d| = d 
will be bigger than the distance do between th e cen ters of 
the two bare harmonic traps, defined in Eq. ( |A2cj ). The 
correction Sx = d — do is the solution of the equation 



d fl o 2 A 

= — I -uuj x + -. —. 

dx I 2 \x + do\ 



/j,u> 2 x 



A 



[x + do) 2 

(A3) 



10 



(it has been assumed do + x > 0) , and can be written as 

(A4) 



Sx = — - sinh 2 1 — In 
3 I 6 



V + 1 + VvW+2) 



with r] = \/[2fiuj 2 (do/3) 3 ]. Expanding to first order in 
e, we find 5x « edo/2. Taking into account that, if the 
traps are sufficiently far apart, the relevant coordinate 
range is X2 > Xi, one obtains the multipole expansion 



A 



|r + d| 



1 - 



ly 2 



lz 2 



d d 2 2d 2 2d 2 

P k (x,y,z) 



E 

k=3 



d k 



(A5) 



where each of the multipole terms P k is a polynomial of 
k th degree in x, y and z - for instance, 



Pa(x, y, z) = -x 



Pi(x,y,z) 



x 2 -°-(y 2 + z 2 ) 



(A6a) 



3x 2 {y 2 + z 2 ) + ^{y 2 + z 2 ) 2 . (A6b) 



By virtue of Eq. ( |A3| ), the linear term in the expansion 
Eq. ( |A5| ) cancels exactly with the one arisin g from the 
harmonic potential /iw 2 (r + d — do) 2 /2 in Eq. (Alt). We 
define the unperturbed Hamiltonians 



p2 i 

H° R = + -Mlo 2 K 2 , H° = H x + H ± , 

R 2M 2 



H r 



where 

Px i 1 2 2 
P i +Pl , 1 2 / 2 , 2 A 

the highcr-multipoles contribution 

A ^ P k {x,y,z) 



E 



<i fe 



(A7) 

(A8a) 
(A8b) 

(A9) 



the force terms 



F(t) = ^(ft 1 +fl a )^(t), /W = -^^m 

a w 2 

(A10) 

where II" is the projector onto the internal state a of 
particle i: 



U x = 11)^11(8)12, 



n 2 = li <8> |l>o <1| ; 



the rescaled frequencies 
v = uj^/I + e, 



= ^V 1 e / 2 ; 

and the shifted coordinate and energy scales 

A 1 



A = R - (1,0,0), E ='- + -^ 2 6x 2 



(All) 
(A12) 

(A13) 



Shifting the coordinate system by Rq — d and the energy 
scale by Eq, we finally obtain Eqs. (13a) and (|13b|). 



2. Fidelity 

The goal of this Section is to show the derivation of 
the analytical temperature dependence of the fidelity, 
Eq. (|3^). We begin by writing the two-particle phase 
as 

tff - yf)) = <W [(-lA + 1] e + o[(a*/d)% 

(A14) 



where 



(Al5a) 



We define 



W § (t) i El +E2 ~ 2 V / - E i^2Cos[w(i 1 - * a 



2fc s r}. 



(A15b) 



S(a,6) = |(XcilXciW)| 2 

= 2(a{l - cos[(k - l)e}} - 6cos[(k - l)e]^A16) 
x (a + b - 1) + (6 - l) 2 + b [b + 2a cos(2ke)] , 



I 2 , b = M 2 , 
|Xci), in the form 1 = J2 a ,p \ c ap\' 
account. From Eq. ( |3l| ) it follows 



where a = |coi| 2 , b = |cio| 2 , and the normalization of 

has been taken into 



min((S(a,6))) = (( minS(a,6))) 

{a,b\ \\ {a,b\ II 



(A17) 



which is a constrained minimization problem, with con- 
straints 0<a<l,0<6<l. The solution cannot 
be found by simply equating the partial derivatives of 
S(a, b) to zero, since - as it will be seen at the end of 
this Section - the minimum turns out to be located at 
the border of the region of allowed parameters. There- 
fore we must take a closer look at the problem to find the 
analytical solution. To this end, we evaluate 



9 a S(a,6) = -2e 2 ( K -l) : 



V K — 1 



d b E(a,b) = -2e 2 ( K + l) 2 (a^-± + ±--b) + 

\ k + 1 2 / 



o[e j, 
(Al8a) 
o[e% 
(A18b) 



While looking for the minimum, we will neglect o[e 4 ] oc 
(a^/d) 8 ^ then, to evaluate it , we w ill use the exact form 
of 2(a, b). According to Eq. ( Alga ), in the region of the 
parameter plane defined by the condition 



1U-1 

b< a- - , 

2 J k+ 1' 



(A19) 



it is db^.{a,b) < — 2k(k +1) < 0. Since k > 1, the 
inequality Eq. (A19) also implies b < 1/2. Therefore the 
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mini mum must be found outside the region defined by 
Eq. (|A19|) , i.e. for 

K+ 1 1 , . s 

a < b + -. (A20) 

K — 1 2 

The latter, by Eq. (|Al8b|) , implies <9 Q S(a,6) < 0. Sum- 
ming up, the minimum is reached for the values (ao,6o) 
of the parameters, where 

o = l, ^ ^— t < b < 1. (A21) 

Z ft + 1 

The problem is therefore reduced to a one-dimensional 
constrained minimization: We have to study the equation 

= 9fcS(ao,6) = 2(cos(2ke) — cos[(k — l)e] 

- 26{cos[(k + l)e] - 1}) , (A22) 

which has solution 

l)e/2] 



- _ 1 sin[(3«; 



(A23) 



2 sin[(rc + l)e/2] ' 

Since the constraint bo < 1 has to be fulfilled, we obtain 
bo = min{6, 1}. Indeed, it is b < 1 only for e > f K (e), 
where 



f K (e) = 2arcsin 
Hence, to o^a^/d) 4 ] 



sin(3Ke/2) — 2shi(ke) 
a/5 + 4cos(2ke) 



(A24) 




«cos[(3K-l)e]»} 
4((cos(Ke)[cos(e) — cos(k£) 



otherwise. 

(A25) 

Now, we f' K (0) = fc/3 > 1 (the inequality following from 
a-oj <C e?). Hence for small T - such that ((e)) <C 1, i.e. 
fceT <C fkj(d/a LU ) 2 -, it is e < / K (e) and the form o f F c i i s 
given by the second row at right-hand side of Eq. ( A25), 
which we can expand in Taylor series around T = for 
taking the thermal average to finally obtain Eq. ((3l|). 
When k — 0, i.e. the third-order anharmonic correction 
is suppressed, 

3(o, b) =\ + 2{a + b- l)(a + b)[l- code)}. (A26) 

The function to be minimized depends now only on the 
sum a + b. Therefore the minimum can be searched by 
fixing one of the two parameters and varying only the 
other one. We can choose ao = 1 as before, and minimize 
Eq. ( A26| ) with respect to b. Eq. (A25) gives the solution 
also in this case. In particular, since f K =o(e) = < e Ve, 
the analytical expression for the fidelity is no w given by 
the first row at right-hand side of Eq. ( A25 ), which for 
k — becomes simply 

1 



F' 



«=o 



: [l + «cos(e) 



(A27) 



The same procedure can be used for the minimization 
over the possible in tern al states also in the q uant um case, 
as is done in App. B 6. By expanding Eq. ( A27) to low- 
est nonzero order in powers of a^/d, we obtain finally 
Eq. 



APPENDIX B: QUANTUM CALCULATION 

In this Appendix we compute both analytically and 
numerically the evolution of the two-ion system, evaluate 
the resulting phase shifts and derive an accurate expres- 
sion for the fidelity, giving as well the explicit expression 
of some quantities used in the text. 



1. Undoing single-particle phases 

In this Section we show how to get rid of the spuri- 
ous phases accumulated during gate operation, in order 
to be left with the gate phase Eq.(||). In the ideal case 
where the external degrees of freedom factorize out at the 
end of the computation, the evolution operator Eq. (0) 
induces both two-particle conditional phases and single- 
particle kinetic phases, depending on each ion's external 
state. When both particles are in their external ground 
state, we can undo the kinetic phases and other inessen- 
tial phases through single-bit operations of the form 



Si = 



E 



a)j(a\ e o , 



(Bl) 



where <pj has to be equal to the kinetic phase acquired 
after the gate operation by particle j, and we choose 



Under these conditions, the compound operator 



U{t) = {S 1 ®S 2 )U{t,t ) 



(B2) 



(B3) 



implements the transformation Eq. (|l|), with d given by 
Eq. (|). 



2. Unperturbed solution 

As explained in the text, since we assume (i > a u we 
can treat the higher-multipoles term H\ as a small per- 
turbation with respect to the rest of the Hamiltonian. In 
this Section we solve exactly the unperturbed problem, 
i.e. we calculate the time-dependent evolution dictated 
by H{t) — Hi. So we want to solve the Schrodinger equa- 
tion 



= [H R (t) + H° r (t) - f(f) ■ (r + d)] !*(*)), 



with initial condition 



(B4) 



(B5) 
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The subscript R (r) denote s the center-of-mass (relative) 
motion, as defined in App. A 1 . The solution is mn 



\Mt))n = e-^ t - t °>/V*< F(t ' )dt ' /S (B6a) 



x exp 



dsK R (s,t ) 



dK* R {s,t ) 



ds 



\Mt)) r 



xe -iK R (t,t )&l e -iKl(t,t )a R \ <t j }R (t Q )) R , 
e -iX?(t-t )/h e idJt f{t')dt'/h ^ mh ^ 

x exp — / dsK r {s,to) 



to 



xe 



-iK r (t,t )al -iK*(t,t )a r 



ds 

\A(to)) r , 



where d R (a r ) is the annihilation operator for the x com- 
ponent of the center-of-mass (relative) motion, and 

K R (t,t Q ) = J_ f dt'F^a^'-^, (B7a) 
V2Mnuj Jt 



K r (t,t ) 



1 



v2fMu J to 



dt' f{t')e iu{t '- to \ (B7b) 



Now the explicit fo rm of the fo rce t erm, Eq. (|22|), can b e 
inserted into Eqs. ( |B7a| ) and flB7b| ) through Eq. ( |A13| ), 
to yield 

i^(*,t ) - (ni + n 2 )jr(w,i,to), (B8a) 
tf r (t,t ) = (ft 2 -ni)«>>Mo)> (B8b) 



where 



K(u,t,t ) = z£-urrSe-W>'l(w,lf) 
l(w,t) = Erf f -- — 



We take as initial state 



(B9a) 
(B9b) 



|^(t„)) fl = |N) fl = |A^) X |JV y ) y |A Z ) Z , (BlOa) 



|Vv(*o)) r = |n) r = 1^)^ |n y > (r^ 



(BlOb) 



During time evolution, the two ions will acquire a state- 
dependent phase shift 



(*(t)|*(to))^|(*(i)|*(M)l^ Q 

which turns out to be given by 



k<*f3 _ 



t>f + <t>f, 



(Bll) 



(B12a) 



« (a + /?) 2 £ 2 $H + (a + /3)VW£XoK, 
-(iVx+iVr +iVzM*-*o), ( B12 b) 
^ pb (a-/3) 2 e 2 $(^)-(a-/3)V^^d/a w 

- [n x v + (n y + n z )u±] (t - t ) , (B12c) 
A 2 



$(w) 



Wr— D(t') 



to 



(B12d) 



with 



.4 
C 



9f[J (a;,* )]-iJ* (w,0), 
«»[J(w,td)], 



2,M 2 



e (WT/2 J 



(B13) 



The equality in Eqs. pl2rj ) and ( |B12c| ) is ap prox imate 
since the integrals in the exponent of Eqs. ( B6a ) and 
(B6b) have been evaluated by means of a saddle-point 



approximation, giving a very good agreement (relative 
difference less than 10 -5 with typical parameters as used 
here) with the exact result, which cann ot be evalu ated 
analytically. Finally, from Eqs. (B6a) and (B6b) we 
obtain 

0{"*>(f,*o)= |(*(*)|tf(* ))| =0 { R rH} 0^\ (B14) 
where 

{ R ,H} = M(-N x aAKR(t,t )\ 2 )e-i\ K ^\ 2 ,(B15 a ) 
- M(-n x , l,\K r (t, *o)| 2 )e-* |jK ' ( * , * 0>|a , (B15b) 
and M{a, b, z) is the confluent hypergeometric function. 

3. First-order perturbation theory 

Now we want to evaluate the lowest-order corrections 
that appear when the higher multipole contributions in 
the Hamiltonian are taken into account. Following pi , 
we expand the evolution operator as 



(B16) 



where Uq (*, to) is the operator o f the unperturbed evolu- 
tion, already calculated in App. B2, and 



Uj (*, to) 



1 



dt \ 



j to J to 



dt 



j-r ■ ■ 



dhU (t,t 3 ) 



x H x U (tj , tj-i) H x U„ (*,_! , tj_ 2 ) (B17) 
x...U Q (t 2 ,t 1 )H 1 U (t 1 ,t ). 

We are interested in evaluating the diagonal matrix ele- 
ments _(\&_(*o)| U(t,t n ) |^(*o)) to first order, acc ordin g to 
Eq. ( |B16| ). Since (U Q (t,t )) is given by Eq. ( pljj) , we 
just need to compute 

(Ui(t,to)) = \ t dt'Offrt'^e^W+f^V^ 
in .h„ 



o i<t>^{t,t ) ft 



ih 



dt'Of(t,t',t ) 



(B18) 
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where the u nperturbed phase factorizes since (as shown 
in Sec. V A ) it does not depend on the initial state, and 
we have defined 



Of(t,t',t ) = K*(io)|E^t,Offi£^(f,to)l*(to)>l 



(B19) 



The exact result, given by Eqs. (B6a) and (B6b), cannot 
be integrated analytically over time. Instead we adopt 
the adiabatic approximation, i.e. we assume that the 
condition (i) of Sect. |y] is satisfied. The Hamiltonian 
then changes slowly enough so that the system, being 
in a motional eigenstate at t = to, follows the changes 
being in the corresponding eigenstate at every subsequent 
time t. This means in our case that, if to < 0, t > 
and their absolute values are large enough, we will have 
\^(t)) « \^>(to)). The relative-motion wavefunction of 
the evolved state is then 



KVr(*o)|Co(t,f)|r>| 



|(r|£/o(*',io)|Vv(*o)>| (B20) 

VVc (a- J "(t) / '^ 2 )^n y {y)^nA z ), 



where e.g. ipn x { x ) = (x\n x ) x G IR. Finally we obtain 



(tfi(t,to))«*e*? fl (A Q/3 + A') 



where 

A af3 = 

A' = 
6 k = 

Si = 



d 

fc=3 



t-tnX 



fe=3 



(B21) 

(B22a) 
(B22b) 



dt'(n\ 



v / 7"~(a^)' c J to 
(n\P k (x,y,z)\n) r 



V2£ 



a^v (l-e)3/4' 



Pk(x + -^a^T(t),y,. 
-P k (x,y,z)\\n) r , (B22c) 
(B22d) 

a v = y/hjjw. (B22e) 



From Eqs. (B16) and (B21) it follows that, to first order, 



(U(t,t )) 



(l7 (*,*o) + ^i(*,*o)> 
(U (t,t )) 1 



\(Uo(t,t ))\ 



, (B23) 



which is equivalent to Eq. @, given that |A Q/3 + A'| = 
\{Ui(t,t ))\ < \{U (t,t ))\ w 1. The internal-state- 
independent part A' cancels out when computing the 
gate phase Eq. (||), as well as the terms of odd k in A a/3 , 
due to the summation over the internal states. The adi- 
mensional quantities 5k and S' k do not depend either on 
the internal state nor on time, but just on the relative 
motional state. We will now calculate them for k = 3,4. 
To be precise, we should not use the eigenstates |n) r of 



H®, as is done in Eq. (B22d), but rather those of the 
full Hamiltonian H r . However, as we will demonstrate 
in the next Section, the corrections are of o[(a w /d) 3 ] and 
therefore we will consistently not take them into account 
in the present calculation. The relevant matrix elements 
are 

(n\x\nf) = -j= (5 n >,n-iVn + 6n',n+iVn~+ l) , (B24a) 

2 r 



ft 

~2 

+ 6 n >, n+2 y/(n + l)(n + 2) 

^3 



On' ,n— 2 \/n (n - 1) + 5 n /,„ (2ra + 1) 

(B24b) 



(n|ir 3 |?i') = 5n',n-3\Ai (« -!)(«■- 2) 

+ 3^, n _in + 3^ )B+ i (n + 1) 3/2 (B24c) 



>/(n + l) (n + 2) (n + 3) 



(n|*V) = ^{j n -, n _4A/n (n-l)(n-2) (n - 3) 

+25 n%n _2(2n + 1)V« («- 1) 
+35 n >, n [2n(n + 1) + 1] (B24d) 
+25 n / jn+2 (2n + 3)V(n+l)(n + 2) 

V(n + 1) (n + 2) (n + 3) (n + 4)}, 



Hence 

4 = 

6' 3 = 
<5 4 = 



1 3 

7= - —5 [2ns + 1 - + n z + 1)] , (B25a) 

vo 2^ 

1 3 

- + [2n x + 1 - v(n y + n z + 1)] , (B25b) 

1 v2^ 



(B25c) 



- [2n x (n x + 1) + 1] - -v(2n x + \){n y + n z + l) 
3 

+ y^ 2 [n y (3n y + 5) + n z (n z + 5) + 4(1 + n y n z )] 

(B25d) 



where v = v /v±_ . 



4. Perturbative corrections to the eigenstates 

Since in our case the perturbation H\ is static, its 
effect on the initial eigenstates of the system must be 
taken into account. In this Section we show how to do 
that in second-order perturbation theory. Our problem is 
to compute the eigenstates of the initial relative-motion 
Hamiltonian 

H r (t ) = H° r +Hi = H° r + eHl , (B26) 

whereby the external force is vanishing at the initial time, 
and 



H, 



l _ hw(f_^ P k (x,y,z) 

w fc=3 



(B27) 
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Therefore we make a perturbative expansion in the small 
parameter e. So we write the eigenstates of H r (omitting 
throughout this Section the subscript r as 



|n(e)) = £V|n«), 



(B28) 



where the first terms are 



|„( 2 )> 



E 



(m(°)|£LV°>) 



F (o) _ R (o) 



(B29a) 



v ( m y^^)( i( yi°? | m (o)- 



l,m^n 



-tin — J (-frn — £■] 



(oh /" 17(0) 



(B29b) 



and |n(°)) are the eigenstates of H®, with eigenenergies 
The k th term in ff r x gives a contribution of order 
~ (hu}/2)(auj/d) k ~ 2 . Since we want to neglect corrections 
of order o[(a w /d) 3 ], we need to go up to k = 4 in the 
expansion of H}. But from Eqs. ( A6a ) and ( A6b ) it 
is straightforward to see that {m^\P$ t ±{x, y, z)\n [u> ) for 
m 7^ n. It follows that 



\n(e)) = \n^)+o[(a u /d)% 



(B30) 



and therefore, as already anticipated in the previous 
Section, for the purpose of the present calculation we 
can consistently use the eigenstates of the unperturbed 
Hamiltonian H2- 



5. Numerical computation 

The goal of this Section is to transform the Schrodinger 
equation for the two-particle wavefunction into a system 
of first-order differential equations for the time depen- 
dence of its projections over the initial eigenstates, better 
suitable for numerical handling. Since the problem has 
cylindrical symmetry around the x axis, the transverse 
coordinates always appear as powers of p = \Jy 2 + z 2 . 
Thus the original three-dimensional problem is equivalent 
to a two-dimensional one. We expand the wavefunction 
(omitting for simplicity the subscript r) as 



OO , r , t 

\m) = E ^nl(t)exp\- d / f(t')dt' 
nJ=0 1 11 L Jt ° 



h[nv + lv± + l)t 



\nl) ,(B31) 



where \nl) = \n) x \l) ± , the \n) x (\l) ± ) are the eigenstates 
of H x (H ± ). From Eq. (|b|) it follows 



Ti 



x (nl\[f(t)x-H x }\n'V) 



^=f(t) (V^e^c^, + V^+Te- wt c n+lil ) 



(B32) 



where the coefficients Cl, correspond to the k term in 



Eq. (A9) - in particular, 

C$ = ^n(n-l)(n-2)e a » t c n ^i (B33a) 
+3n 3 / 2 e wt c„_ 1 , i + 3(n + l) 3 / 2 e^ t c n+hl 
+ y/(n + l){n + 2)(n + 3)e~ l3,yt c n+3 j 

+ ^/n(l + l)(l + 2)e^- 2 ^ t c n ^ u+2 
+ (21 + 1) y^e ivt c n - x , x + V^TTe- wt c n+u ) 
+ v/(n + 1)1(1 - \)e-^- 2v ^c n ^ 2 
+ y/(n + l)(l + + 2)e~^ +2 ' J ^ t c n+u+2 



C$ = 



-y/n(n - l)(n - 2)(n - 3)e 4wt c„_ 4 j (B33b) 
+ y/n(n - \)e 2wt [3D [^/iQ^e 2 ^^,^ 
+ ^/(l + l)(l + 2)e- 2 ^ t c n ^ l+2 
+ {21 + l)c„_ 2 , ; ] - 2(2n - l)c„_2,l} 

if) 2 



■^l(l-l)(l-2)(l-3)e 4 ^ t c n , l - 4 



+Z0y/l(l-\) 



(2n + l)--(2i-l) 



2ivxt„ 
e C„;_2 



9£> 2 

— [21(1 + 1) + 1] - 3z>(2n + 1)(2Z + 1) 



3 . 



+3 [2n(n + 1) + 1] |c ni - -uy/Q + 1)(Z + 2) 

x [z/(2Z + 3) - 4(2n + l)]e^ 2 ^*c n!i+2 
- ^ + W + m + m + 4)e- 4 ^ t c„, i+4 
+ v /(n + l)(n + 2)e- 2wt {3 I >[(2Z + l)c n+2il 
+ ^/(l + l)(l + 2)e- 2i ^ t c n+2 , l+2 
+ y/l(l - l)e 2 ^ t c„ +2 , i _ 2 ] - 2(2n + 3)c n+2 ,z} 
-y/(n + l)(n + 2)(n + 3)(n + 4) e - 4i,yt c„ + 4j. 

Excitations higher than a certain level should be ab- 
sent as long as we are in an adiabatic regime. Thus in 
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Eq. ( B32j ) we neglect the coefficients above a certain N. 
We have checked that the result is independent of the 
cutoff. 



6. Fidelity 

The goal of this Section is to evaluate the gate op- 
eration fidelity in the full three-dimensional quantum- 
mechanical framework. The overall phas e B ap p earing 
in E q. (^5| ) can be computed from Eqs. ( B12a )-( B12d ) 
and ( |B22a| )- p22d| ), as 



where the minimization over the co effic ients {c Q ^} has 
been carried out exactly as in App. [A 2| . Here, only the 
relative motion comes into play because Ad is indepen- 
dent of the center-of-mass motion, and cos( Ad) has been 
expanded up to o(A8 3 ). Hence Eq. ( phf ) follows, by tak- 
ing into account that e <C 1 and 8 C \ f=a 8 = tt, and ex- 
panding in a Taylor series for hu; <C kgT. 



6 w 2u{ v^£t 



X 



o OO , 



k=3 



{t - to)|, 
(B34) 



where 5'^ is defined in Eq. ( B25d| ) , and it has been taken 
into account that e <C 1- In the ideal case, according 
to Eq. (|l|) for -8 = tt, the gate operation transforms the 
initial internal state |x) into 



7. Many-ions calculation 



In this Section we simply give the definitions of the 
parameters appearing in Eq. (p0[): 



1x0 = E (-l) Q/5 ^|a) 1 |/?) 2 . (B35) 

a, 0=0 

In a more realistic situation the initial total density op- 
erator <jt at a temperature T is given by 



<j T = p T (t ) ® |x)(x|, 



(B36) 



where pr(to) is defined in Eq. (p6[), and we recall that 
lu w v w After the gate operation we have 

»r= E c c e /34'/3'G^p T (io)G^ /3 ,|a)i(a'|®|/3) 2 (/3'|, 

(B37) 

where G a p = (aP\G\a[5) , and the gate opera tor G is 
defined in Eq. (Esl). As already stated in Sect. V A, be- 



cause of adiabaticity, the motional state after the gate 
operation is unchanged, i.e. G a p PT^o)G^ a ,^, ~ pt(^o)- 
If 6 = it, the minimum fidelity F(T), given by Eq. (^|), is 



6 oo 



F(T) = min(l- 7 ) 6 nE 



7 '{n 

i—1 n.i — 1 
2 i |„ \2\f\„ |2 i \ n |2 



(IcooP + lciil 2 )^ 



2(|coo| 2 + |cii| 2 )(|coi| 2 + |cio| 2 ) cos(A0) 

(hi 



,i| 2 + |c 10 | 2 ) 2 



(l- 7 ) 3 



II E 7" I (« 4 |[l + cos(A0)]|n l )(B38) 



— i 2 - - d 2 H 



Wi = wyl + e?/i 



2 1 + e% ' 



(B39a) 
(B39b) 



1 w 



1 - &. 



i [i/> (2) («) + ^ (2) (^ + 1 - i) + C(3)] ,(B39c) 

1 X -S I 



2 ^ \i - i 
i [^(iV + i-i)-^ 1 ^) 



(B39d) 



where 7ife is the harmonic number and tp^(z) the 
polygamma function of order fc, and £(s) is the Riemann 
zeta function. It is 



max |^| = C(3) « 1.2, max M = — « 0.82. (B40) 

2,n z,n lz 



3/)2 



6 3 



(l + e 5)(i_ e 2/ 4) 



f-r - 

V d J i 



-hu/k B T 



(1 _ e -Hv/k B T\2 



